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ABSTRACT 

Many blazars detected by the Fermi satellite, observed spectroscopically in the optical, 
are line-less, and have been classified as BL Lac objects. Optical-UV photometry of nearly 
one hundred of them allowed to determine the redshift for a handful of objects and redshift 
upper limits for the great majority. A few of these are candidates to be "blue quasars", namely 
flat spectrum radio quasars whose broad emission lines are hidden by an overwhelming syn- 
chrotron emission peaking in the UV. This implies that the emitting electrons have high en- 
ergies. In turn, this requires relatively weak radiative cooling, a condition that can be met if 
the main radiative dissipation of the jet power occurs outside the broad line region. We con- 
firm this hypothesis by studying and modelling the spectral energy distributions of the 4 "blue 
quasars" recently discovered. Furthermore, we discuss the distribution of Fermi blazars in the 
7-ray spectral index - 7-ray luminosity plane, and argue that "blue quasars" objects are a 
minority within the blazar populations. 
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1 INTRODUCTION 

The Fermi satellite is detecting 7-ray emission from a large number 
of blazars (Flat Spectrum Radio Quasars, FSRQs, and BL Lacs). 
From the data collected in the first two years of operation a "clean" 
sample was constructed from the data of the Fermi/Large Area 
Telescope (LAT) instrument, the 2LAC sample (Ackermann et al. 
201 1), that allows population studies. 

The 2LAC sample includes 395 sources classified as "BL 
Lacs", 310 FSRQs, 157 sources of "unknown type", 4 Narrow Line 
Seyfert 1 (Abdo et al. 2009) and other 18 "non blazar" AGNs and 

2 starburst galaxies. Of the 395 sources classified as BL Lacs, 56% 
lack a redshift determination, which limits the possibility of dis- 
cussing their physical properties. When an emission lines is visible, 
the subdivision between the BL Lac and FSRQs categories is based 
on the equivalent width (EW) of the line, as measured in the rest 
frame: the blazar is classified as BL Lac if the rest frame EW of 
any permitted line is smaller than 5 A (Stickel et al. 1991). 

In order to gain redshift information for BL Lacs without any 
visible line, Rau et al. (2012, hereafter R 12) set up a program to ob- 
tain simultaneous photometry over a wide wavelength range using 
the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) 
and the Swi/z/Optical Ultraviolet Telescope (UVOT). 

80 blazars with optical/radio identification but without redshift 
information were selected from the 2LAC sample based on celestial 
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position and small foreground reddening; 8 more (2 with known 
redshift) were included though not part of the clean 2LAC sample 
because of Fermi data quality problems; other 16 2LAC sources 
with known redshift were included for verification. In total, 104 
blazars were considered: 82 have been classified as BL Lacs by 
Ackermann et al. (201 1), 3 as FSRQs and 19 as of unknown type. 

All the sources (but one, due to a lack of precise coordinates) 
were observed simultaneously with GROND and with the Swift 
UVOT. GROND can observe simultaneously in 7 filters, from 2.4 
fim (K band) to 4000 A (g band), while UVOT can observe in 6 
filters (one filter at the time), from 5400 A (v band) to ~ 1600 A 
(uvw2 filter; with center wavelength at 2000 A). These data, cover- 
ing the wavelength interval from 2.4 /im to 1600 A, allow to derive 
photometric redshifts if the source is far enough to be affected by 
intervening Lya absorption. 

As a result, photometric redshifts were determined at the 90% 
confidence level for 1 1 sources, including 3 blazars that also had a 
previously measured spectroscopic redshift (in one case it is only a 
lower limit on z, derived by the presence of an intervening absorp- 
tion line). Furthermore, the absence of any Lya absorption feature 
in the spectrum could lead to the estimate of an upper limit on the 
redshift for 81 blazars, including 12 blazars with known spectro- 
scopic redshift. The upper limits in these 12 sources were all con- 
sistent with the spectroscopic measurement. 

Among the 11 blazars with z estimated photometrically, 
Padovani, Giommi & Rau (2012, P12 hereafter) discussed 4 
sources, chosen because of the flat (in vF v ) optical continuum 
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(i.e. F v oc v~ a with a ^ 1). They constructed their spectral en- 
ergy distributions (SEDs), including the X-ray data of the X-Ray 
Telescope (XRT) onboard the Swift satellite and the 7-ray data of 
Fermi. They concluded that, despite their featureless optical spec- 
trum and high synchrotron peak frequency, more typical of low 
power BL Lac objects (if blazars obey the "blazar sequence" as 
proposed by Fossati et al. 1998), these sources are probably FS- 
RQs whose broad emission lines are swamped by the beamed non- 
thermal continuum. 

The aim of the present paper is twofold. The first aim is to an- 
alyze in detail the SEDs of the 4 blazars discussed in PI 2 in order 
examine whether the large 7-ray luminosity together with a high 
frequency synchrotron peak is in contrast with the physical inter- 
pretation of the blazar sequence given in Ghisellini et al. (1998). 

This interpretation assumes that the peak frequency of the 
blazar synchrotron emission depends mainly on the energy of the 
emitting electrons which, in turn, depends on the rate of radiative 
cooling. In this scheme, we expect that high power blazars (i.e. 
FSRQs) emit most of their radiation within the broad line region 
(BLR), which can be the dominant source of seed photons for scat- 
tering. Radiative cooling, dominated by the Inverse Compton pro- 
cess on this external population of seed photons (i.e. seed photons 
produced externally to the jet, so the inverse Compon process is 
called external Compton, EC for short) is severe, and therefore the 
energy of the relevant electrons, emitting at the peak of the spec- 
tral energy distribution (SED), is small. At the other extreme of 
the blazar sequence we have low power BL Lacs, that lack (or have 
weak) broad emission lines. EC is much less important, and the cor- 
responding cooling is unimportant. Electrons can then reach high 
energies, and produce high synchrotron (and self Compton, SSC) 
frequencies. 

However, as described in Ghisellini & Tavecchio (2008, see 
also Georganopoulos et al. 2001), intermediate situations may exist 
if the dissipative region of the jet lies beyond the BLR. In this case 
the relativistic electrons are not subject to strong EC losses, and 
the main emission processes become synchrotron and SSC, with 
some contribution from EC scattering of IR photons produced by 
the torus. The resulting SED is then foreseen to be similar to the 
SED of classical BL Lac objects. It is then possible to have a "blue " 
quasar i.e. an object with emission lines and at the same time a 
SSC-dominated jet. These conditions should occur in sources with 
a relatively small accretion disc luminosity id and/or a relatively 
large black hole mass M. 

The second aim of this paper is to take advantage of the upper 
limits in redshift, made possible by the combination of GROND 
and UVOT observations, to locate these blazars in the a 7 -L 7 
plane, to see if they violate the general trend observed in Ghisellini, 
Maraschi & Tavecchio (2009). In that paper, we analyzed only the 
blazars with spectroscopic redshift detected in the first 3 months of 
Fermi. We can now update the sample using the 2LAC catalog and 
including blazars with redshift information provided by R12. 

We use a flat cosmology with Ho = 70 km s _1 Mpc -1 , 
f2M=0.3 and the notation Q — W x Qx in cgs units. 



2 ANALYSIS OF SWIFT DATA 

P12 presented the data of the 4 blazars in their sample only in the 
form of SED plots. Since the Swift data are publicly available on 
the Swift archive, we have analyzed both the XRT and UVOT data 
of the in order to check the results of P12. 

The data were screened, cleaned and analysed with the soft- 



ware package HEASOFT v. 6.12, with the calibration database up- 
dated to 22 March 2012. The XRT data were processed with the 
standard procedure (XRTP IPELINE v. 0.12. 6). All sources 
were observed in photon counting (PC) mode and grade 0-12 (sin- 
gle to quadruple pixel) were selected. The channels with energies 
below 0.3 keV and above 10 keV were excluded from the fit and 
the spectra were rebinned in energy so to have at least 20-30 counts 
per bin in order to apply the \ 2 test - When there are no suffi- 
cient counts, we applied the likelihood statistic as reported by Cash 
(1979). Each spectrum was analysed through XSPEC v. 12.7.1 with 
an absorbed power law model with a fixed Galactic column density 
as measured by Kalberla et al. (2005). The computed errors repre- 
sent the 90% confidence interval on the spectral parameters. The 
X-ray spectra displayed in the SED have been properly rebinned to 
ensure the best visualization. 

UVOT (Roming et al. 2005) source counts were extracted 
from a circular region 5"-sized centred on the source position, 
while the background was extracted from an annulus with internal 
radius of 7" and variable outer radius depending on the nearest con- 
taminating source. Data were integrated with the uvot imsum task 
and then analysed by using the uvot source task. The observed 
magnitudes have been dereddened according to the formulae by 
Cardelli et al. (1989) and converted into fluxes by using standard 
formulae and zero points. 

The Swift data analysed by us are the same used in R12 and 
PI 2. We found substantial agreement in three cases, but for the 
blazar RX J0035.2+1515 we found UVOT fluxes rather different 
from R12, and for this source we found no break. Note that in the 
field of RX J0035.2+1515 there is a very bright star at a distance 
of 36" (TYC 1187-1355-1, with B=10.39) that causes problems 
when subtracting the background. For its estimation, we have used 
an annular region of size 7"-20" centered on the blazar (see also 
§4.1 below). 



3 BL LACS OR FLAT SPECTRUM RADIO QUASARS? 

The 4 considered blazars have a featureless optical continuum, and 
can be classified as BL Lacs if the classical definition is adopted 
(i.e. an equivalent width of the lines less than 5 A). However, P12 
suggested that these blazars are instead FSRQs, whose emission 
lines are swamped by the relativistically boosted jet flux. The main 
argument for this classification is the strong radio power of these 
sources, that is typical of FSQRs. We agree with this interpretation, 
and we would like to offer another argument in favor of the FSRQ 
classification of these blazars. 

Ghisellini et al. (201 1) and Sbarrato et al. (2012), considering 
Fermi detected blazars, have shown that there is a correlation be- 
tween the 7-ray luminosity and the luminosity of the broad lines, 
that includes not only FSRQs, but also sources that are classified as 
BL Lacs according to the classical subdivision (based on the EW 
of the lines). Furthermore, if the luminosities are measured in Ed- 
dington units, there is a divide between BL Lacs and FSRQs for 
LsLn/L-Edd ~ 5 x 10~ 4 and for L 7 /L E dd ~ 0.1. The latter 
values is obtained using the isotropic equivalent of the 7-ray lumi- 
nosity, i.e. the K-corrected 7-ray flux multiplied by (A-ivd^), where 
dh is the luminosity distance. This of course does not imply that L 7 
is isotropic. If the photometric redshift of our 4 blazars is correct, 
then their 7-ray luminosity is large, exceeding the L 7 /I/Edd ~ 0.1 
value even for black hole masses equal to M = 1O 9 M . An ex- 
ception could be RX J0035.2+1515, if the true redshift is z ~ 0.3, 
since in this case L 7 ~ 10 45 erg s _1 , making L 7 /Z/Edd < 0.1 
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for M > 8 x 10 7 M o . In this case we should consider the source 
as a BL Lac, namely a blazar whose emission lines, if present, are 
intrinsically weak. 

The correlation found in Sbarrato et al. (2012) concerns 
sources closely aligned with the line of sight, and it is foreseen 
that when the 7-ray sensitivity will improve, what we see now as 
a correlation is in fact a boundary in the Lblr-L 7 plane. Bearing 
this in mind, the correlation has the form 

L BL r~4L°' 93 (1) 

with a large scatter, since the 7-ray luminosity is highly variable 
in single objects even when averaging over one or two years (see 
e.g. Ghirlanda et al. 2011). This offers a rough way to estimate 
the luminosity of the broad lines. When a good optical spectrum 
is available, we can then suggest the minimum ratio - between the 
boosted non-thermal and the thermal continua - needed to hide the 
lines. 



(Condon et al. 1998) and vL v ~ 2 x 10 42 erg s _1 . The source is de- 
tected in the infrared by the WISE satellite in all its four filters (3.4, 
4.6, 12 and 22 (imfl The corresponding data points are consistent 
with the extrapolation of the spectrum derived from the GROND 
fluxes. 

We have re-analyzed the UVOT data, finding a very bright 
source at ~ 36" from RX J0035. 2+1515, as mentioned in §2. Es- 
timating the background in a region of the sky free of sources, we 
have derived de-reddened fluxes significantly smaller than the ones 
reported in R12, and a harder spectrum, with no sign of a break. 
There is then the possibility that the derived photometric redshift 
is affected by the uncertainties caused by incorrectly subtracting 
the background. For this reason, we will consider for this source 
both the photometric redshift derived by Rll, and also 2 = 0.3. 
This roughly corresponds to the lower limit on z due to the non- 
detection of the host galaxy both in the image and in its possible 
contribution to the SED (see Wagner et al. 1996; Sbarufatti, Falomo 
& Treves 2005). 



4 NOTES ON INDIVIDUAL SOURCES 
4.1 RX J0035.2+1515 

The source has been observed by the SDSS to have a featureless 
continuum. NED reports z — 1.09. On the other hand SDSS reports 
z — 1.057 as a result of an automatic analysis, and also alerts 
that z is actually unknown. The quoted values are not believable. 
Adopting the photometric redshift given in R12, z = 1.28, we 
can derive a (5a) upper limit on the flux of the Mgll line which is 
the most prominent broad line observable in the spectral range of 
SDSS. We derive L Mg n < 1.4 x 10 43 erg s _1 . We then use the 
template given in Francis et al. (1991), adding the Ho contribution 
(not included in Francis et al. 1991), with a relative weight of 77 (on 
a scale in which the Lya is 100). The total weight of all lines is then 
555 (see Celotti, Padovani & Ghisellini et al. 1997), and the weight 
of Mgll is 34. Therefore we derive Lblr. = (555/34) LMgii < 
3 x 10 44 erg s -1 . With a covering factor of 0.1, the upper limit on 
the accretion disc luminosity is Ld < 3 x 10 45 erg s _1 . The 0.1 
value for the covering factor is uncertain and should be taken as 
an average value with some dispersion (see e.g. Baldwin & Netzer 
1978; Smith 1981). 

Assuming a standard, geometrically thin optically thick disc 
(Shakura & Sunyaev 1973), the peak frequency of its spectrum has 
VpL d (v p ) < 1.5 x 10 45 erg s _1 , a factor ~ 20 below the ob- 
served vL v in the optical, that has a luminosity ~ 3 x 10 46 erg 
s . We do not have any telling observation for the value of the 
black hole mass. Sbarrato et al. (2012), studying a sample of FS- 
RQs detected by Fermi, calculated the average black hole mass of 
FSRQs: (M) = 5x 1O 8 M0. Each black hole mass was estimated 
by Shen et al. (2011) through virial methods. Using this average 
black hole mass, the ratio Ld/Lsdd < 0.05. 

As long as id ~ 10 -2 LEdd, the hypothesis of a standard disc 
is justified (radiatively inefficient disc should corresponds to Ed- 
dington ratios smaller than 0.01). 

As discussed in §3, the correlation between the BLR luminos- 
ity and the 7-ray luminosity in the FermifLAT energy band offers 
a way to estimate Lblr and Ld. Using Eq.Q]and L 7 ~ 4 x 10 46 
erg s _1 , we obtain Lblr ~ 10 44 erg s _1 (and Ld ~ 10 45 erg s -1 ) 
with an uncertainty of at least a factor 4. Reassuringly, this estimate 
is consistent with the value found above. 

The radio information are poor, since only the 1.4 GHz NRAO 
VLA Sky Survey (NVSS) point is available, with a flux of 1 8.7 mjy 



4.2 SUMMS J053748-571828 

Not observed by SDSS, its photometric redshift is z = 1.55. To 
estimate the presence or not of a standard accretion disc, and there- 
fore its BL Lac or FSRQ nature, we can use the correlation between 
L 7 and Lblr, giving Lblr ~ 6.3 x 10 43 erg s _1 and there- 
fore an accretion disc 10 times more powerful. Please note that 
the dispersion around the L 7 -Lblr is large, so the above values 
should be taken as an order of magnitude estimate. Nevertheless, 
since the optical continuum in this source has a luminosity simi- 
lar to RX J0035. 2+1515, it is conceivable that the synchrotron flux 
has swamped the (indeed present) broad emission lines, if the op- 
tical spectrum has a S/N similar to RX 10035.2+1515 (i.e. ~ 40), 
or worse. The line of reasoning is the following: assuming that the 
photometric redshift is correct, we know at what frequencies the 
prominent broad emission line (i.e. Mg II) should appear. In order 
to be visible, this line should have a minimum luminosity, depend- 
ing on the quality of the spectrum. Since R12 state that the optical 
spectrum is featureless, we can then assign a lower limit on the line 
luminosity assuming a S/N ratio. We find that a BLR of luminosity 
~ 6 x 10 43 erg is consistent with the absence of lines in the optical 
spectrum if the S/N <40. 

The source is detected by WISE. Although not simultaneous 
with GROND, the IR data points lie on the extrapolation of the 
spectrum defined by the GROND points. 

In the radio, the source flux is 99.8 mjy at 843 MHz, as re- 
ported by the Sydney University Molonglo Sky Survey (SUMSS; 
Mauch et al. 2003). This corresponds to vL v ~ 10 43 erg s" 1 , 
slightly larger than RX J0035.2+1515 at 1.4 GHz. 

4.3 CRATES J0630-2406 

Not observed by SDSS, its photometric redshift is z = 1.6. The 
L 7 -Lblr relation gives Lblr ~ 5.5 x 10 44 erg s _1 , suggesting 
a rather luminous disc (Ld ~ 5.5 x 10 45 erg s _1 , within on order 
of magnitude). If true, this disc luminosity would correspond to 
a VpLd(v p ) a factor 20 below the optical flux: a spectrum of the 
same quality of the SDSS spectrum of RX J0035.2+15 15 would not 
reveal any line. The source is detected by WISE. The corresponding 

1 Cutri et al. 2012: http://wise2.ipac.caltech.edu/docs/release/allsky/ 
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Figure 1. Energy densities as seen in the comoving frame of a function of 
distance from the black hole, -Rdiss- The adopted mass is M = 5xlO 8 M 
and the disc luminosity is L d = 3x 10 45 erg s —1 . These are the values used 
to model the SED of RX J0035+1515 when assuming 2=1.28. The vertical 
grey line marks the value of Rdiss use d in Fig. [5] The jet is assumed to 
accelerate with T oc (-Rdiss /Rs) 1 ^ 2 untu it reaches its maximum value 
(in this case V = 18). At -Rdiss ~ l- 7 X 10 17 cm, the source is outside 
the BLR, and at -Rdiss ~ 3 X 10 17 cm is of the same order of U' g 

and f/j'p, if the IR radiation produced by a molecular torus intercepts 10% 
of Ld- Further out, but within Rip, the radiation energy produced by the 
torus dominates again the radiative cooling, and the spectrum is dominated 
again by the EC process. The label "jump" indicates the jump in external 
radiation energy density as we move outside the BLR. 



IR data well match the extrapolation of the spectrum defined by the 
GROND points. 

The large radio luminosity, exceeding 10 44 erg s _1 , indicates 
a powerful source, suggesting that this blazar is a FSRQ. 



4.4 CRATES 1312-2156 

This source can be identified with PKS 1309-216. The absorption 
line observed by Blades, Murdoch & Hunstead (1980), identified as 
CIV, gives a lower limit on the source redshift: z > 1.491. The ra- 
dio luminosity is large, exceeding 10 44 erg s^ 1 . This suggests that 
this source is a FSRQ, whose accretion disc continuum and lines 
are swamped by a strong synchrotron component. We can estimate 
(approximately) an upper limit to the contribution of the disc, by 
assuming that, if observed with the SDSS, broad lines do not ap- 
pear. With a quality of the spectrum similar to the SDSS spectrum 
of RX J0035.2+15 15, and assuming the photometric redshift, the 
most prominent observable broad line would be Mgll. Since it is 
not observed, the accretion disc must lie well below the observed 
continuum, with a peak VpL^ifp) at least a factor ~ 20 below the 
optical continuum, i.e. below 5 x 10 45 erg s _1 . This value is of the 
same order of what derived by the L 7 -Lblr correlation, that gives 
L d ~ 6 x 10 45 erg s _1 . Also this source is detected by WISE. 
Although not simultaneous with GROND, the WISE data lie on the 
extrapolation of the spectrum defined by the GROND points. 
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Figure 2. Top panel: The entire SED of 0035+1515, assuming 2 p j 1 =1.28. 
The (green) empty circles are the UVOT data as analyzed by us. The black 
line in the optical is the (de-reddened) SDSS spectrum. The WISE not si- 
multaneous data (labelled "WISE not sim") are also shown. The solid line 
is a synchrotron plus EC model, the thick black dashed line is a pure SSC 
model. The grey dotted line is the contribution of the assumed accretion 
disc, IR torus and X-ray corona, as labeled. The dashed line is the con- 
tribution of the SSC flux for the SSC+EC model. Bottom panel: The SED 
assuming z = 0.3. The shown model is a pure SSC model, assuming no 
external radiation. 



5 MODELING 

The adopted model is a one-zone and leptonic model, fully de- 
scribed in Ghisellini & Tavecchio (2009). The main properties are 
summarized in the Appendix, mainly to explain the meaning of the 
parameters listed in Tab. Q] and Tab. [2] In extreme synthesis, the 
model assumes that the bulk of the jet dissipation takes place in 
one zone located at some distance i?diss from the black hole. For 
simplicity, the emitting region is assumed spherical with a radius 
R = ip -Rdiss, with tp = 0.1. The region is moving with a bulk 
Lorentz factor F, and is observed under a viewing angle 9 V . Ener- 
getic electrons are injected throughout the source for a time equal 
to the light crossing time R/c, and the particle distribution is cal- 
culated (though the continuity equation) at this time, considering 
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Figure 3. Top panel: SED of RX J0537-5718 (top); CRATES 0630-2406 
(mid) and 1312-216 (bottom). The model assumes the presence of external 
radiation, as in Fig. [2] For CRATES 0630-2406 we also show a pure SSC 
model (thick black dashed line). 



radiative losses and possible electron-positron pair production and 
their reprocessing. 

A specific feature of our modeling is that it takes into account 
how the relevant energy densities change as a function of iidiss, the 
distance from the black hole of the main dissipation region. Fig. [T] 
illustrates this point, showing how the different contributions to the 
energy density (as seen in the comoving frame) change as a func- 
tion of i?diss (see the lower axis), and, equivalently, as a function of 
^diss/^s (top axis), where Rs is the Schwarzschild radius (for a 
black hole of mass M — 5 x 1O 8 M0). Since the jet is assumed to 
accelerate (with T = min[Tf, (i?diss/3i?s) 1//2 ]), the bulk Lorentz 
factor is relatively small in the inner regions of the jet. As a conse- 
quence, the external radiation is not strongly boosted. On the other 
hand, if the Poynting flux initially dominates the total power of the 
jet, the inner regions will be magnetically dominated. As a result 
the main source of cooling is synchrotron and SSC. When F be- 
comes large, but still 7?diss < -Rblr, the main source of cooling 
is the inverse Compton scattering off the broad line photons (seen 
strongly enhanced in the comoving frame). The adopted scaling be- 
tween i?BLR and the disc luminosity implies that, in the comoving 
frame, the external radiation energy density simply is: 



4tt j R| lrC 



12tt 



erg cm 



(2) 



where we have assumed that Lblr = O.lLd and -Rblr = 
10 17 .£^ 45 cm. Beyond Rblr this energy density decreases fast, 
because it is not enhanced any longer by the Doppler boosting. 
There is, however, another important external component, if the 
source has a molecular torus re-emitting in the infrared the in- 
tercepted disc radiation. Again, we assume a scaling between the 
size of the torus, Rm, and the disc luminosity of the kind Rm = 
2.5 x 10 18 Ly, 2 r , cm. Assuming that the torus intercepts a fraction 
a of the disc radiation, as long as i?di ss < ^?ir the corresponding 
energy density measured in the comoving frame is: 
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where a_i = a/10 - . Comparing Eq.[2]with Eq. [5]we see that, 
going from inside to outside the BLR there is a "jump" of two or- 
ders of magnitudes in external radiation energy density, as shown 
in Fig.[U as i?diss moves from just inside to just outside -Rblr- 

Note that the magnetic energy density becomes comparable 
with t/gLR and Urn only in a narrow range of distances, just be- 
yond Rblr- 

The different cooling rate reflects in a different electron en- 
ergy distribution AT (7), having a break at a different 7 C defined as 
the random Lorentz factor of those electrons loosing half of their 
energy in one dynamical time: 



7c = 



3m c c 



(4) 



Aa T RU' 

If these are the electrons responsible for the peak of the SED, we 
expect a large difference between blazars emitting inside or outside 
the BLR. In fact, the jump in U' is only partly compensated by the 
increase of R. Assume a change of a factor 3 in J?diss (between be- 
ing inside and being outside the BLR), and a corresponding change 
of a factor 3 in the magnetic field (B oc R^ 1 ). Since the jump in 
U' is nearly 100, we expect that the peak synchrotron frequency v s 
of the region R ou t is larger than the R[ n one by a factor: 



^s,out -^out7c,out ^ou 

V S ,in -Bin7 2 in B in 



RinU[ n 



Ro 



300 



(5) 
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where the subscripts "out" and "in" stand for outside and inside 
the BLR. It follows that if a "standard" FSRQ with R diss < Rblr 
has a synchrotron peak at u s ~ 10 Hz, then a blue FSRQ - with 
-Rdiss just a factor 3 larger (sufficient to make -Rdiss > Rblr) - has 
a synchrotron peak at ~ 3 x 10 15 Hz, in the UV. The example above 
is only illustrative, since the different cooling will affect the shape 
of the particle distribution, not only its break, but it does explain 
the expected difference between "red" and "blue" FSRQs. 

5.1 Physical parameters of blue FSRQs 

The results of the models are shown in Fig. [2] and Fig. [3] and the 
model parameters are listed in Tab. [Tj and [2] We alert the reader 
that since our model assume a one-zone emitting region, that must 
be quite compact, we cannot account for the radio emission, that 
must be produced in much larger regions, not to be self-absorbed. 
The observed self-absorbed synchrotron frequency of our model is 
typically ~1000 GHz. Another concern is the value of the black 
hole mass. The model could estimate it if the accretion disc con- 
tribute to the optical-UV flux, but in all cases we found instead that 
it is dominated by the non-thermal beamed continuum. Therefore 
we simply assumed a black hole mass value similar to the average 
value (i.e. M~5x 10 s M Q ) found in Sbarrato et al. (2012). 

RX J0035.2+1515 and CRATES J0630-2406 — We have applied 
a SSC+EC model and a pure SSC model for both sources. In the 
first case the model assumes the presence of a standard, Shakura- 
Sunyaev disc and a reprocessing torus intercepting 10% of the disc 
radiation, and reemitting it in the IR. The accretion disc luminos- 
ity has been chosen to agree with the limits given in the previous 
section. The radiation produced by the torus, together with the (not 
greatly boosted) BLR emission is the dominant component for the 
inverse Compton scattering producing the very high energy flux, 
while the SSC component gives a comparable contribution in the 
X-ray band. In these cases -Rdiss > Rblr, but only by a fac- 
tor 1.5-2.5. The total (magnetic plus radiative) energy density as 
seen in the comoving frame is small, and the energy of the elec- 
trons emitting at the peaks of the SED is large. The synchrotron 
peak frequency is beyond 10 15 Hz for both sources, and the syn- 
chrotron flux overwhelms the disc emission. The magnetic field is 
not far from equipartition with the radiation energy density, imply- 
ing roughly equal synchrotron and Self Compton luminosities. 

We have then applied a pure SSC model, to see if the presence 
of the disc and torus is required or not. As can be seen, this model 
can account for the IR-7-ray SED quite satisfactorily, at the cost 
of increasing the Doppler factor S (achieved by increasing F and 
decreasing 8 V ), and lowering the magnetic field. 

Finally, we have assumed 2 = 0.3 for RX J0035.2+1515 and 
we have applied a pure SSC model, to see if with this redshift the 
source could be fitted with the same parameters of high frequency 
BL Lacs. Indeed, we find consistent values, very similar to other 
sources of this kind (see Tavecchio et al. 2010; Ghisellini et al. 
2010). From the point of view of model fitting, we cannot decide if 
RX J0035.2+15 15 is a blue FSRQs or a high frequency peaked BL 
Lacs, but in any case we require a large 7 pca k and a small comoving 
energy density (see the discussion below). 

SUMMS J05 3748-57 1828 and CRATES 1312-2156 — For these 
two sources we confirm the break in the bluest UVOT filters. In 
the absence of absorption, this break would be inconsistent with 
the extrapolation of the X-ray spectrum at lower frequencies. We 
have applied only one model, assuming the presence of a standard 
accretion disc and torus, and z — z p h- With this redshift, the radio 
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Figure 4. The Lorentz factor of the electrons emitting at the synchrotron 
peak, 7 P eak, as a function of the radiation plus magnetic energy density, 
as seen in the comoving frame. The 5 sources discussed in this paper are 
labelled and are compared to a collection of blazar samples, discussed in 
Celotti & Ghisellini (2008) (CG); in Ghisellini et al. (2010, 2011) and in 
Sbarrato et al. (2012) (G10+G1 1+S12). As can be seen, the "blue quasars" 
studied here lie along the main trend of decreasing 7 pca k when increasing 
the energy density. 



luminosity is large: this suggests that these sources are powerful 
FSRQs. Again, we have assumed a disc luminosity in agreement 
with the limits/indications found in the previous section, a torus 
re-emitting 5-10% of Ld, and Rdiss nearly twice -Rblr- The bulk 
Lorentz factor (T — 12) and viewing angle (9 V — 3°) are rather 
usual for Fermi FSRQs (Ghisellini et al. 2010). The assumption of 
Rdiss > Rblr implies a large 7 pca k and a small comoving energy 
density. 

Concerning the jet power (see Tab.|2j, our sources have values 
similar to the other Fermi blazars (studied in Ghisellini et al. 2010), 
but certainly not extreme, at the border between BL Lacs and FS- 
RQs. Consider that since the average electron energy is large, there 
is little difference in considering the presence of one (cold) proton 
per emitting electron or not. Even a pure electron-positron pair jet 
would carry a similar total jet power. 

The main finding of the modeling is that, in order to obtain 
a good agreement with the data, a large 7 pea k and a small energy 
density t/ B + V' x is required in any case, either if the model is a 
pure SSC or if it includes external radiation. Fig. |4]compares the 
obtained values of 7 pca k and f/ B + U[ with the ones of other 7-ray 
detected blazars. As can be seen, all our sources follow the same 
trend of other blazars. This implies that they are consistent with the 
explanation of the blazar sequence based on radiative cooling, as 
suggested in Ghisellini et al. (1998). 

Rau et al. (2012) considered about 100 blazars, finding an up- 
per limit on the redshift for the majority of them, and a photometric 
value for 11. Of the latter, only 4 (the ones considered by PI 2) 
have a flat optical shape (spectral index a < 1, F(y) oc ^~ Q ). 
The remaining blazars are then either closer or steeper. Therefore 
blue quasars (powerful blazars with a < 1) are rare objects, as pre- 
dicted in Ghisellini & Tavecchio (2008) on the basis of the small 
parameter space occupied by these sources. Consider also that if 
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Name 


2 


Pkliss 


M 


PtBLR. 


P{ 


L d 


B 


r 


6>v 


7b 


7max 




S2 


Tpcak 


U' 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


m 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


0035 (EC) 


1.28 


300 (2000) 


5e8 


173 


5.0e-4 


3.0 (0.04) 


1.28 


18 


1.5 


7e3 


2e5 


-1 


3 


7561 


0.1137 


0035 (SSC) 


1.28 


300 (2000) 


5e8 




7.0e-4 




0.18 


24 


1.3 


1.5e4 


3e5 


1.6 


2.8 


1.18e4 


2.1e-3 


0035 (SSC) 


0.3 


180 (1200) 


5e8 




4.0e-4 




0.08 


20 


2 


3.2e4 


5e5 


1.5 


4 


3.8e4 


4.4e^t 


0537-5718 


1.55 


390 (2600) 


5e8 


229 


3.2e-3 


5.25 (0.07) 


1.39 


12 


3 


6e3 


1.8e5 


-0.5 


2.7 


5583 


0.114 


0630-2406 


1.60 


375 (2500) 


5e8 


194 


5.2e-3 


3.75 (0.05) 


0.86 


17 


2.2 


Le4 


2.5e5 





2.7 


1.16e4 


0.0553 


0630 (SSC) 


1.60 


312(1300) 


8e8 




1.7e-3 




0.38 


22 


1.5 


3e4 


2.5e5 


0.75 


3. 


2.07e4 


0.0102 


1312-216 


1.77 


480 (2000) 


8e8 


245 


0.014 


6 (0.05) 


0.84 


14 


3 


3e3 


le5 


-1 


2.2 


4484 


0.0601 



Table 1. List of parameters used to construct the theoretical SED. Not all of them are "input parameters" for the model, because Pblr. is uniquely determined 
from Ld. Col. [1]: name; Col. [2]: redshift; Col. [3]: dissipation radius in units of 10 cm and (in parenthesis) in units of Schwarzschild radii; Col. [4]: black 
hole mass in solar masses; Col. [5]: size of the BLR in units of 10 15 cm; Col. [6]: power injected in the blob calculated in the comoving frame, in units of 
10 45 erg s — 1 ; Col. [7]: accretion disc luminosity in units of 10 45 erg s — 1 and (in parenthesis) in units of LEdd; Col. [8]: magnetic field in Gauss; Col. [9]: 
bulk Lorentz factor at Pdiss; Col. [10]: viewing angle in degrees; Col. [11] and [12]: break and maximum random Lorentz factors of the injected electrons; 
Col. [13]: and [14]: slopes of the injected electron distribution [(3(7)] below and above 7^; Col. [15]: random Lorentz factors of the electrons radiating at the 
synchrotron peak; Col. [16]: sum of the radiation and magnetic energy density in the comoving frame. The total X-ray corona luminosity is assumed to be in 
the range 10-30 per cent of L^. Its spectral shape is assumed to be always oc u~ 1 exp(— hu/150 keV). 



Name 


logP r 


logP B 


logP c 


logPp 


0035 (2=1.28, EC) 


44.20 


45.25 


43.09 


43.66 


0035 (2=1.28, SSC) 


44.05 


43.80 


44.26 


45.72 


0035 (2=0.3, SSC) 


42.91 


42.49 


43.61 


44.59 


0537-5718 


44.65 


45.20 


43.48 


44.12 


0630-2406 


45.16 


45.05 


44.12 


44.56 


0630 (SSC) 


44.85 


44.41 


44.28 


44.39 


1312-216 


45.42 


45.08 


44.42 


45.02 



Table 2. Logarithm of the jet power in the form of radiation (P r ), Poynting 
flux (Pb), bulk motion of electrons (P c ) and protons (P p , assuming one 
proton per emitting electron). Powers are in erg s — 1 . 

RX J0035.2+1515 is indeed at more moderate redshifts, it could be- 
come a high frequency peaked BL Lac, implying that blue quasars 
are even rarer. In turn this implies, in the framework of our model, 
that the dissipation region in the jet is almost always within the 
BLR. 



6 THE UPDATED FERMI BLAZARS ' SEQUENCE 

Ghisellini, Maraschi & Tavecchio (2009) pointed out the existence 
of a well defined trend in the « 7 -L 7 plane of the blazars de- 
tected after the first three months of Fermi observations. BL Lacs 
and FSRQs were separated in this plane, with BL Lacs having the 
flattest slopes. These were very bright blazars, detected at more 
than the 10a level. Now the 2 years catalog is available, and we 
can compare the locations of the R12 sources in the updated ver- 
sion of the a 7 -L 7 plane. Fig. [5] shows the results considering all 
BL Lacs (blue empty points) and all FSRQs (filled red circles) 
with redshift present in the clean 2LAC sample. What plotted is 
the [0.1-10 GeV] K-corrected luminosity. Overplotted to these 
sources, we show the blazars considered by R12 that have spec- 
troscopic and photometric redshifts, together with the blazars for 
which the photometric analysis yields only a upper limit on the 
redshift (2 < 2 max )- For the latter sources, we have considered a 
range of redshifts, between 0.3 and 2 max , hence a range of L 7 . The 
choice of 0.3 as the minimum possible redshift follows from the 



fact that, at smaller redshifts, the host galaxy starts to be visible. 
Note the following: 

• Comparing this figure with the analogous figure in Ghisellini, 
Maraschi & Tavecchio (2009), the trend between a 7 and L 7 is now 
even more clear. Likely, this is due to the improved software anal- 
ysis, assigning to the background more low energy photons than 
before. As a result, steep spectra (a 7 > 1) become flatter, while 
already flat sources are less affected. We recall that the observed 
a 7 -L 7 trend is partly due to selection effects: when deeper surveys 
will be available, we will start to see slightly misaligned blazars of 
lower L 7 (and same « 7 ) as well as blazars with smaller black hole 
masses (and presumably with smaller L 7 ), so that the trend will be 
much more dispersed. 

• R12 consider some blazars with spectroscopic redshift, but 
with no redshift indication in the 2LAC catalog. These are plotted 
with different symbols in Fig. [5] There are two groups: five blazars 
have a flat a 7 and are candidates to be high frequency peaked BL 
Lacs. The remaining ones are located in the region populated by 
FSRQs, and are then candidates to be normal FSRQs. 

• The blazars with the estimated photometric redshift have all 
large L 7 . The 4 blazars studied here are the flattest. Remember that 
they were selected to have a flat optical spectrum: they also have a 
flat 7-ray spectrum. Fig.[5]shows how the blazar RX J0035.2+1515 
would be located if its redshift were z — 0.3. In this case it would 
be a normal highly peaked BL Lac object. 

• The blazars in R12 with only a limit on z are located in the re- 
gion of the plot already populated by sources with known 2. There 
are no outliers. The existence of very powerful BL Lac objects with 
flat 7-ray spectra is excluded. This agrees with the blazar sequence. 

We can conclude that the observed boundary in the a 7 -L 7 
plane at large luminosities is not due to selection effects, namely 
to the difficulty to determine the redshifts of powerful BL Lac ob- 
jects. It must have a physical origin. We (Ghisellini, Maraschi & 
Tavecchio 2009) have proposed that this is due to blazars with the 
largest black hole masses, accreting at different rates. 

Small accretion rates correspond to weak discs: the ionizing 

luminosity L lon is small, the broad line region is correspondingly 

1/2 

weak and very close to the black hole (if -Rblr oc L ion still holds). 
Dissipation takes place beyond the BLR, the radiative cooling is 
weak, electrons reach high energies, and from an observational 
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Figure 5. The energy spectral index a 7 as measured by Fermi as a function of their [0. 1-10 GeV] (rest frame) luminosity. Filled (red in the electronic version) 
points are FSRQs, empty (blue) circles are BL Lacs. The sources studied in R12 are plotted with larger symbols if the redshift is determined (spectroscopically 
or photometrically). Sources with only un upper limit for the photometric redshifts are plotted assuming a luminosity the upper limit as the maximum redshift, 
and z m i n = 0.3 as the minimum (this corresponds approximately to the limit at which the host galaxy should be visible). For RX J0035. 2+1515 we plot also 
the location of the source assuming z = 0.3. 



point of view we have a high frequency peaked BL Lac with no or 
very weak broad lines (being the synchrotron peak in the UV/X-ray, 
it is easier for the jet emission to overwhelm the weak disc and BLR 
emission). Increasing the accretion rate, we enter in the regime of 
radiatively efficient discs and then a normal BLR. If the jet dissi- 
pates mainly within the BLR, the radiative cooling is strong, imply- 
ing a low frequency peaked FSRQs. Furthermore, if the jet power is 
proportional to the accretion rate M, then L 7 is also linked with it, 
and therefore we have the chain: low M — > low jet power + weak 
and small BLR — > flat a 7 + small L 7 . And the opposite for large 



M. This chain can be repeated for each black hole mass. Of course, 
the lowest the black hole mass, the lowest the maximum L 1 . 



7 CONCLUSIONS 

We have studied the nature of a few blazars, whose redshift has 
been recently determined, photometrically, by Rau et al. (2012). In 
addition, we have considered all blazars for which Rau et al. (2012) 
could establish a lower limit on the redshift. Our conclusions agree 
with Padovani et al. (2012): the 4 blazars considered in our and their 
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papers are probably FSRQs (with a residual doubt about the deter- 
mination of the photometric redshift of RX J0035. 2+1515, because 
the UVOT data analysis is complex). 

Independently on the actual redshift of RX J0035. 2+1515, all 
4 blazars are "blue" (i.e. the synchrotron and the inverse Compton 
peak frequencies are large) because the radiative cooling is weak. If 
they do have broad lines, they are overwhelmed by the strong syn- 
chrotron emission. This agrees completely with the standard ex- 
planation of the blazar sequence. The fact that it was possible to 
estimate a photometric redshift for a few sources, out of a sample 
of nearly one hundred, confirms that these blue FSRQs are rare. 
This in turn implies that it is rare, for a jet, to produce most of the 
radiation we see at large distances from its black hole. 

The found blue quasars can have accretion discs emitting at a 
level of a few per cent of the Eddington one, "normal" black hole 
masses (namely between 10 8 and 10 9 Mq) and bulk Lorentz fac- 
tors. Their properties are similar to the average properties of other 
blazars. 

We have considered the updated q 7 -L 7 plane, considering all 
sources with spectroscopic redshifts within the 2LAC sample and 
including the blazars for which Rau et al. (2012) could estimate an 
upper limit on z. We find no outliers in this plane: in other words, 
the upper limits on the redshift exclude the presence of very pow- 
erful BL Lacs. This confirms our previous explanation of the sepa- 
ration of BL Lacs and FSRQs in this plane based on a sequence of 
accretion rates. 
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APPENDIX 

At a distance Pdiss from the black hole of mass M the jet dissipates 
part of its power and injects relativistic electrons throughout the 
emitting region, assumed to be spherical, with radius R = ^ Pdiss, 
with ip — 0.1. In the region there is a tangled magnetic field B. The 
relativistic electrons are injected with a smoothly joining broken 
power law in energy: 



Q(7) = Qo 



1 + (7/7b)~ sl+s 



(6) 



The energy particle distribution ^(7) [cm -3 ] is calculated solving 
the continuity equation where particle injection, radiative cooling 
and pair production (via the 7-7 — > e process), are taken into 
account. The created pairs contribute to the emission. 

The injection process lasts for a light crossing time R/c, and 
we calculate N(-y) at this time. This assumption comes from the 
fact that even if injection lasted longer, adiabatic losses caused by 
the expansion of the source (which is traveling while emitting) and 
the corresponding decrease of the magnetic field would make the 



observed flux to decrease. Therefore the calculated spectra corre- 
spond to the maximum of a flaring episode. 

The total power injected into the source in the form of rel- 
ativistic electrons is P( — m c c 2 V J Q(j)jdj, where V = 
(4/7r/3)i? 3 is the volume of the emitting region. 

The bolometric luminosity of the accretion disc is Ld ■ Above 
and below the accretion disc, in its inner parts, there is an X-ray 
emitting corona of luminosity Lx (it is fixed at a level of 30% 
of Ld). Its spectrum is a power law of energy index ax = 1 
ending with a exponential cut at E c =150 keV. The specific en- 
ergy density (i.e. as a function of frequency) of the disc and the 
corona are calculated in the comoving frame of the emitting blob, 
and used to properly calculate the resulting External inverse Comp- 
ton spectrum. The BLR is assumed to be a thin spherical shell, of 
radius Pblr = 10 17 L d ^ 2 5 cm. We consider also the presence of a 
IR torus, at larger distances. The internally produced synchrotron 
emission is used to calculate the synchrotron self Compton (SSC) 
flux. Table[T]lists the adopted parameters. 

The power carried by the jet can be in the form of radiation 
(P r ), magnetic field (Pb), emitting electrons (P e , no cold electron 
component is assumed) and cold protons (P p , assuming one proton 
per emitting electron). All the powers are calculated as 



P = irR 2 r 2 l3cUi 



(7) 



where (/.■ is the energy density of the i component, as measured in 
the comoving frame. 

The power carried in the form of the produced radiation, P r = 
7rP 2 r 2 /3c?7 r ' ad , can be re-written as [using (7 r ' ad = L' / (4/7rP 2 c)]: 

r 1 



r 2 r 2 

4 4<5 4 



45 2 



(8) 



where L is the total observed non-thermal luminosity (1/ is in the 
comoving frame) and I7 r ' ad is the radiation energy density produced 
by the jet (i.e. excluding the external components). The last equality 
assumes 9 V ~ 1/r. 

When calculating P c (the jet power in bulk motion of emitting 
electrons) we include their average energy, i.e. U' c = n c ("f)m c c 2 . 
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